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20.  Abstract  (Continued) 


Progress  for  the  period,  1  September  1977  to  31  August  1978,  includes:  (1)  development  of  a 
method  for  calculating  energy  deposition  profiles  from  range  distributions,  (2)  development  of  a 
method  to  estimate  energy  deposition  in  ion  bombardments  which  incorporate  transport  of  energy 
downstream  by  primary  knockon  atoms,  (3)  an  investigation  of  the  influence  of  1  atomic  percent 
Al,  Mo,  Si  and  Ti,  respectively,  to  the  swelling  behavior  of  binary  nickel  alloys  under  neutron 
irradiation;  (4)  an  investigation  of  swelling  in  refractory  metals  and  alloys  (Mo,  Nb,  V,  TZM,  Nb-1% 
Zr  and  V-20%  Ti)  at  625°C  for  fast  neutron  fluences  to  5.4  x  1022  n/cm2;  (5)  an  investigation  of 
deuterium  trapping  in  radiation  damaged  aluminum  and  stainless  steel  using  elastic  scattering  tech¬ 
niques  to  measure  the  increase  in  dislocation  density  produced  by  radiation  damage  in  copper  single 
crystals. 
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Program  Description 

The  Cooperative  Radiation  Effects  Simulation  Program  (CORES)  was 
initiated  voluntarily  by  five  Branches  from  the  Engineering  Materials  (now  the 
Material  Science  and  Technology  Division)  and  the  Radiation  Technology 
Division  of  NRL  on  the  basis  of  their  common  interests  in  the  problems  of 
simulating  radiation  damage  in  metals.  The  program  promotes  the  exchange 
of  information,  discussion  of  problems,  and  the  pursuit  of  collaborative 
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interest  to  the  damage  simulation  problem.  The  major  portion  of  the  work  is 
sponsored  by  the  Office  of  Naval  Research.  Since  research  findings  which 
apply  to  the  objectives  of  one  sponsor  may  also  be  of  interest  to  others,  the 
overall  progress  related  to  damage  simulation  is  included  in  the  written  report. 
Several  of  the  participating  Branches  have  independent  programs  on  other 
aspects  of  the  radiation  damage  problem;  when  results  obtained  in  these 
programs  are  judged  to  be  of  interest  to  CORES  participants  they  may  also  be 
included,  informally,  in  the  CORES  program  review. 
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COOPERATIVE  RADIATION  EFFECTS  SIMULATION  PROGRAM 


Annual  Progress  Report 
for  the  Period 

1  September  1977  -  21  August  1978 


SUMMARY 

I.  RADIATION  DAMAGE  THEORY 

A.  Use  of  Range  Distribution  to  Approximate  Energy  Distributions 

Range  distributions  can  be  used  to  estimate  the  distribution  in  energy 
which  develops  in  an  ion  beam  as  it  penetrates  a  target.  Several  variants  of  this 
approach  are  discussed  and  applied  to  the  calculation  of  energy  deposition 
profiles. 


B.  Approximation  for  Energy  Deposition  in  Ion  Beam  Bombardment 

A  method  is  developed  for  estimating  energy  deposition  in  ion  bombard¬ 
ment.  By  making  extensive  use  of  a  table  of  range  versus  energy,  it  is  possible 
to  estimate  the  energy  transport  effects  due  to  primary  knock-on  atoms 
depositing  energy  downstream  from  the  sites  of  their  creation. 

II.  HEAVY  ION  AND  NEUTRON  DAMAGE  STUDIES 

A.  The  Effect  of  Alloy  Additions  on  the  Swelling  of  Neutron  Irradiated 
Nickel 

The  void  and  dislocation  structures  in  pure  nickel  and  the  binary  alloys  Ni-l 
at.%  Al,  Ni-l  at.%  Mo,  Ni-l  au9&  Si,  and  Ni-l  at.%  Ti  were  examined  following 
neutron  irradiation  to  2.2  x  l(j  n/cnri  (E  >0.1  MeV)  at  455°C.  The  dislocation 
structures  of  all  the  materials  were  similar,  consisting  of  tangles  of  line 
dislocations,  many  of  which  intersected  voids.  Comparing  the  void  formation  of 
the  alloys  to  that  in  the  pure  nickel,  the  following  observations  were  made:  the 
Al  addition  did  not  significantly  affect  the  void  formation,  except  to  make  it 
more  homogeneous  and  to  narrow  the  size  distribution;  the  Mo  and  Ti  both 
increased  the  swelling  by  about  a  factor  of  two,  with  the  Mo  increasing  the  void 
number  density  and  the  Ti  increasing  the  void  size;  the  Si  decreased  the  swelling 
by  about  a  factor  of  ten,  by  reducing  the  void  size.  These  results  differ 
significantly  from  those  of  an  earlier  ion  irradiation  study  of  the  same  alloys,  in 
which  all  of  these  alloying  elements  reduced  or  eliminated  the  void  formation. 

The  reason  for  these  differences  is  most  likely  related  to  different  void 
nucleation  mechanisms  being  important  in  the  two  types  of  irradiation. 
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B.  The  Microstructure  of  Neutron  Irradiated  Refractory  Metals  and 
Alloys 

The  refractory  metals  Mo,  Nb,  and  V  and  the  alloys  TZM,  Nb-l%Zr  and  V- 
20%Ti  were  examined  by  transmission  elytron  microscopy  following  irradiation 
at  650°C  to  fluency  of  5.4  x  n/cm  ,  E>0.1  MeV  (all  materials  except  V- 
20%Ti)  and  8.4  x  10^  n/cm  ,  E>0. 1  MeV  (Mo  and  V-20%  Ti).  Voids  representing 
significant  volume  swelling  were  observed  in  the  Mo  (2.3  -  3.0%),  TZM  (4.0%),  Nb 
(1.5%)  and  V  (3.6%).  Some  precipitation  occurred  during  irradiation  in  all 
materials  except  the  Mo  and  V.  The  results  indicate  that  void  swelling  can 
probably  be  controlled  at  moderate  levels  at  this  temperature  in  refractory 
alloys.  Degradation  of  mechanical  properties  due  to  precipitation  of  interstitial 
impurities  may,  however,  be  a  more  serious  problem. 

III.  LIGHT  ION  BEAM  STUDIES 

A.  Deuterium  Retention  in  Aluminum  and  Stainless  Steel  With  and 
Without  Helium  Loading 

Thin  samples  of  aluminum  (1.2  mils)  and  stainless  steel  (0.6  mils)  were 
implanted  with  deuterium  both  with  and  without  prior  implantations  of  helium. 
These  samples  were  then  profiled  to  determine  concentrations  of  helium  and 
deuterium  by  elastic  scattering  techniques.  Analysis  of  these  measurements 
show  deuterium  was  retained  in  aluminum  foils  both  with  and  without  previous 
helium  implantations  but  deuterium  was  retained  in  stainless  steel  only  with 
prior  helium  loading  indicating  enhanced  deuterium  retention  with  predamage  in 
stainless  steel. 

B.  X-ray  Diffraction  Measurements  of  Low-Level  Radiation  Damage 

X-ray  diffraction  gives  a  quantitative  measure  of  eprly  stages 
damage  in  metals.  Alpha-particle  irradiations  of  2  x  10  ,  4  x  10  , 

C/crr>  in  popper  jingle  crystals  show  an  apparent  threshold  of  little  damage  up 
to  2  x^lO"^  C/crrr  and  then  a  nearly  linear  increase  of  dislocations  up  to  6  x  10'^ 
C/cm  .  Preliminary  experiments  with  polycrystalline  copper,  not  discussed  in 
this  report,  seem  to  show  similar  results. 


of  radiation 
and  6  x  10"^ 


COOPERATIVE  RADIATION  EFFECTS  SIMULATION  PROGRAM 


Annual  Progress  Report 
for  the  Period 

1  September  1977  -  21  August  1978 


INTRODUCTION 


F.  A.  Smidt,  3r.,  Coordinator, 

Thermostructural  Materials  Branch 
Material  Science  &  Technology  Division 

This  report  summarizes  annual  research  accomplishments  of  the  inter- 
divisional  cooperative  program  now  in  its  eighth  year.  The  goal  of  this 
program  has  been  to  use  ion  simulation  techniques  to  increase  the  understand¬ 
ing  of  neutron  damage  of  materials  for  advanced  nuclear  systems.  Light  and 
heavy  ion  bombardment  techniques  are  used  to  represent  neutron  damage  in 
order  to  accelerate  research  and  to  permit  an  evaluation  of  critical  param¬ 
eters  which  might  not  be  feasible  using  nuclear  reactors  as  test  facilities.  In 
this  program  the  techniques  applied  have  advanced  along  with  our  growing 
understanding  of  radiation  damage  of  materials.  To  complement  the  ion 
damage  simulation  are  other  programs  such  as  theoretical  evaluations  of 
atomic  collisions,  computation  of  energy  deposition,  and  parallel  experiments 
using  nuclear  reactors.  The  productivity  of  the  CORES  Program  is  illustrated 
by  a  full  list  of  technical  papers  and  presentations  in  the  last  section. 


Note:  Manuscript  submitted  July  12,  1979 
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RESEARCH  PROGRESS 


I.  RADIATION  DAMAGE  ANALYSIS 

A.  Use  of  Range  Distributions  to  Approximate  Energy  Distributions 

(Irwin  Manning,  Radiation-Matter  Interaction  Branch, 

Radiation  Technology  Division) 

Introduction 


Consider  a  uniform  ion-beam  of  energy  E„  bombarding  an  amorphous  - 
target  in  slab  geometry,  with  the  beam  parallel  to  the  x-axis  and  the  target 
surface  located  at  x  =  0.  We  are  often  interested  in  the  single-particle  flux 
density  ^(e,x)  which  gives  the  distribution  in  energy  which  develops  in  the  ion 
beam  when  it  penetrates  to  a  distance  x: 

<f>  (E,x)  dE  =  number  of  beam  atoms  with  energy  dE  about  E 

which  cross,  per  unit  bombarding  flux  per  unit  time, 
a  perpendicular  unit  area  at  the  depth  x.  (1) 

According  to  this  definition, 0is  normalized  to  unit  bombarding  flux: 


/ 0(E,O)  dE 


=  1. 


(2) 


Estimates  of  this  flux  density  are  difficult  to  obtain,  but  estimates  are 
much  more  readily  available  for  the  range  distribution  f(x,  Eg), 

f(x)  dx  =  number  of  beam  atoms  per  unit  perpendicular 
area  per  unit  bombarding  flux  per  unit  time 
coming  to  rest  at  the  depth  dx  about  x.  (3) 


For  example,  in  addition  to  this  range  distribution,  reference  1  provides  for  a 
large  variety  of  beam-target  combinations,  tables,  giving  as  a  function  of  beam 
energy,  the  total  range  R  (Eg),  the  projected  range 


R 


<Eb) 


<  X  >  = 


oo 

I 


f(x)  dx, 


(4) 


and  the  range  straggling a(Eg),  where 


2 


00 

a  2  (Eb)  =  <(x  -  x)2  >  =  J (x  -  Rp)2  f(x)  dx. 


(5) 


According  to  the  above  definition,  the  range  distribution  is  normalized  to  unity: 


f(x)dx  =  1. 


(6) 


This  paper  concerns  a  class  of  methods  for  using  the  range  distribution  f(x)  to 
approximate  the  energy  distribution  <f>  (E,  x).  Usually  one  will  utilize  the  esti¬ 
mate  as  the  integrand  inside  an  integral,  which  will  tend  to  smooth  out  errors  in 
*(E,  x).  Section  VI  below  discusses  the  employment  of  this  estimate  to  obtain 
energy  deposition  profiles. 


Energy  Distribution  from  the  Range  Distribution 

We  ignore  the  angular  spread  which  develop  in  the  penetrating  beam  as  it 
comes  to  a  stop,  assuming  that  the  particle  velocities  are  always  parallel  to  the 
x-axis.  Consider  the  flux  of  particles  f(x2»dx2)  which  will  stop  at  dx,  about  X2* 
We  make  a  model  in  which  there  are  no  range^fluctuations  (in  the  ordinary  sense) 
among  beam  particles.  In  reality,  different  beam  particles  undergo  different 
fluctuations  in  range,  but  on  the  average  the  beam  obeys  the  range-stopping 
function  R  (E).  These  fluctuations  are  removed  in  our  model,  and  instead, 
fluctuation?  are  introduced  into  the  model  by  assigning  different  range-stopping 
functions,  R  (E,  xJ  to  different  beam  atoms,  picked  in  such  a  way  that  the  model 
has  the  rangi  distribution  ftx,).  For  the  above  flux  of  atoms  which  will  stop  in  dx2 
about  x,,  consider  their  energy  at  some  x.<  X2*  The  fluctuations  in  our  model 
are  sucn  that  all  of  these  atoms,  when  at  x  =  Xj,  lie  within  a  small  energy 
interval  dE,  about  the  energy  E..  There  is  a  one-to-one  correspondence  between 
those  beam  atoms  which  have  energy  Ej  at  x  equal  to  Xj  and  those  atoms  which 
will  stop  at  x  equal  to  X2*  We  have 

<t>  (Ej,  Xj)  dEj  =  f(x2)  dx2,  (7) 

where  Ej  is  determined  by  Xj  and  X2  (and,  implicitly,  Eg): 

x2  =  x1+RP(E1,x2)  (8) 

which  gives  Ej  as  a  function  of  Xj  and  x2: 

Ej  =  Ej  (Xj,x2).  (9) 

Once  a  specific  model  giving  the  altered  range-stopping  function  R  (E,  X2)  is 
chosen,  E.  (x.,  x,)  is  determined  by  equation  (8)  and  equation  (7)  then  gives  an 
estimate  of  the  energy  distribution.  For  a  well-behaved  function  Ej  (xj,  x2),  this 
estimate  is  exact  in  the  limit  Xj-»  X2,  so  we  expect  this  estimate  to  be  a  good 
one  for  *2~  "small". 
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k  *k  °nC  eas*ly  verif‘es  that  the  above  estimate  (7)  obeys  an  identity  satisfied 
by  the  true  energy  distribution.  Consider  the  particle  current 

jW  =  E,T)dE,  (10) 

where  0  is  the  reactor  flux  corrsponding  to  scalan  flux  (1).  Gauss’  law  expresses 
the  conservation  of  beam  particles  as 


divj  =  -p(r),  (u) 

where  p  is  the  volume  density  of  stopping  particles.  For  our  slab  geometry, 
equations  (10)  and  (11)  become  y’ 


d/<?x  j(x)  =  -f  (x) 


where 


j(Xj) 


y* dE10(x1, 


x) 


(12) 


(13) 


all  allowed  E 


The  estimate  (4)  gives 

00 

Xxj)  =  J' dx2  f(x2), 

o 

which  obeys  the  identity  (12). 


(14) 


We  now  go  on  to  consider  specific  models  for  R  (E,  xj  and, 

\nnmolu  P  »  I  P  7  2  J 


correspondingly,  E^Xj.xJ. 

Case  A:  No  Fluctuations  for  x  >  x 


Here  we  assume  that  different  beam  particles  have  different  range- 
stopping  functions  Rp  (Ej,  x?)  only  for  x<  Xj.  For  x  >  Xj  we  have,  for  all  x2, 


RplEj,^)  =  Rp(Ej), 

Equation  (8)  becomes 


(15) 


X  >  Xj. 


(16) 


As  will  be  discussed  in  Section  VI,  this  expression  is  equivalent  to  what  was 
used  in  references  6  and  4. 
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Case  B:  Uniform  Stretching 
Here  we  choose 

R  (E,  x2)  =  C(x2)Rp(E)  (17) 

with  the  constant  C  chosen  such  that 

x2  =  C  Rp  (Eg)  (18) 

Equation  (8)  than  gives 


VEl’  =  '‘(V^VV  <19> 

as  will  be  discussed  in  Section  VI,  this  expression  is  equivalent  to  what  was 
used  in  references  7,  8,  and  9.  The  corresponding  energy  distribution  was  used 
by  Marwick  (8)  to  calculate  PKA  recoil  spectrum. 


Case  C:  Fluctuations  Proportional  to  Range  Straggling 
Here  we  choose 

R  (E,  x2)  =  R  (E)  +  C(x2)a(E), 
with  the  constant  C  chosen  such  that 
x2  r  Rp  (Eg)  +  Ca  (Eg). 

In  this  case  the  model  fluctuations  are  proportional  to  range  straggling 
R  (E,  x2)  -  R  (E)  =  Ca(E) 

Equation  (8)  becomes 

x2  =  +Rp(E1)  +  ^a(E1)/«(EB))((x2-Rp(EB)). 


(20) 


(21) 

(22) 


(23) 


Energy  Deposition  Profiles 

As  a  beam  particle  traverses  the  target,  it  deposits  energy  into 
elastic  collisions  given  by  (2). 


CO 


Sn(E)  = 


(dE/dR)elaS.ic  -  /Td» n- 


(24) 


where  T  is  the  target  atom  recoil  energy,  and  da  is  the  atomic  section  for 
elastic  collisions.  The  corresponding  quantity  for  (aE/dx)  can  be  approximated 
as 
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(25) 


(dE/dxelastic>  Sm<E)  ^R(E)/dRp). 

The  profile  S  (x)  for  energy  deposited  into  elastic  collisions  by  the  entire  ion 
beam  as  it  traverses  the  depth  x  can  be  written 


Sg^x) 


J dE  g(E,x)  Sn(E). 


(26) 


A  more  refined  approximation  can  be  obtained  by  employing  the  expression 
(25)  (Reference  3): 


sB2(x) 


f  g(E,x)  Sm(E)  (dR(E)/dRp^. 


(27) 


A  still  more  refined  approximation  can  be  obtained  by  including  the  LNST 
theory  for  energy  partition  (5)  to  obtain  (3) 

E„ 


sD(x) 


B 

f  dEg(E,x)  SL(E)fdR(E)/dRpV 

o  '  ' 


where 


SL(x)=  N  f  T  ?7  (t)  d  <j 


(28) 


(29) 


and  T}  (T)  is  the  LNST  partition  function. 


On  using 
becomes: 

the  relation  (7),  the  above  energy 

dispostion  profile 

SB(xi)  = 

OO 

/dx 2  f(x2) 

X1 

OO 

(30) 

sB2(xi)  = 

J  dx2  f(x2)  Sm  (Ej)(dR(E)/dRp\  , 

(31) 

and 

x  '  ' 

X1 

sd(xi>  = 

OO 

f  dx2  f(x2)  SL  (Ej)  (dR(E1)/dRp\ . 

(32) 

X 

XI 

For  case  A,  SB  of  equation  (29)  becomes  the  energy  deposition  profile  calculated 
by  Kulcinski,  Laidler,  and  Doran  (6),  while  Sg-  and  SD  of  equations  (31)  and 
(32)  become  the  quantities  calculated  by  Manning  and  Mueller  with  the  computer 
code  E-DEP-1  (4).  For  case  B,  S~  of  equation  (32)  becomes  the  energy  deposition 
profile  calculated  by  Johnston,  Rosalowski,  Turkalo  and  Lauritzen  (7),  as  described 
by  Marwick  (8)  and  further  disucssed  by  Matthews  (9). 
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Marwick  (8)  and  Matthews  (9)  point  out  that  for  case  A  one  obtains 
a  non-physical  energy  distribution  at  the  target  surface  whereas  for  case  B, 
one  correctly  obtains  E  =  E„  as  the  only  allowable  energy.  These  results  are 
readily  confirmed  by  examination  of  equations  (16)  and  (19)  at  the  point  x  =  0. 
We  see  from  equation  (23)  that  case  C  also  has  the  correct  energy  distribution 
at  the  surface. 
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B.  Approximation  for  Energy  Deposition  in  Ion  Beam  Bombardment 

(Irwin  Manning,  Radiation- Matter  Interaction  Branch, 

Radiation  Technology  Division) 

Introduction 

As  do  some  other  energy  deposition  calculations  for  ion  beam  bom¬ 
bardment,  the  NRL  Code  E-DEP-I  makes  the  approximation  that  a  primary 
knockon  atom  (PKA)  deposits  its  energy  close  to  the  site  of  creation  (4).  For  low 
bombarding  energies  this  approximation  can  lead  to  large  errors,  since  the  PKA  in 
actuality  deposits  energy  downstream  from  its  creation  site.  Brice  has  formulated 
a  method  for  taking  those  recoil  effects  into  account  (10).  The  present  work 
formulates  a  different  approach  for  calculating  these  effects  by  making  extensive 
use,  in  the  spirit  of  Reference  (11),  of  a  table  of  range  distributions,  the  present 
work  makes  different  (and  probably  milder)  approximations  than  does  Brice. 

Notation  for  Elastic  Collisions 


As  in  Reference  (11),  we  consider  a  beam  of  ions  with  atomic  number  Z, 
and  energy  E^  bombarding,  in  slab  geometry,  an  amorphous  target  of  atoms  with 
atomic  number  Z ^  To  distinguish  binary  collisions  between  beam-target  atoms 
from  binary  collisions  between  target-target  atoms,  we  shall  denote  the  former  by 
subscript  12  and  the  latter  by  22. 

We  shall  write  the  elastic  cross  section  as  K^EjT),  where  E  is  the  energy 
of  the  (Z.)  beam  atom,  and  T  is  the  recoil  energy  of  the  (Z^)  target  atom.  We 
have,  for  The  penetrating  beam 


K^EjT)  dT  dR  =  number  of  Z^  recoils  of  energy  dT  about  T  produced 
per  unit  time  as  an  ion  beam  of  unit  fluence  and 
energy  E  traverses  a  path  length  dR  through  a  unit 
perpendicular  area.  (1) 


For  an  ensemble  of  penetrating  ions  Zj,  we  also  have 

Kj^  (E,T)  dTdR  =  ensemble  average  of  the  number  of  recoils  of  energy 
dT  about  T  produced  as  a  penetrating  ion  Zj  of 
energy  E-traverses  a  path  length  dR  (2) 

For  the  elastic  energy  loss  we  have 

S",2(E>  '  <dE/dR)elaStic  *  N / dl  TK12(E,T)  (3) 


where  N  is  the  density  of  target  atoms. 
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Let  0  (E,x)  be  the  energy  distribution  of  the  penetrating  (Z^)  beam 
defined  in  Reference  11  and  let  *  (T,x)  be  the  distribution  function  for  PKA's 
(that  is,  recoils  produced  in  12  collisions)  produced  by  this  beam: 

*  (T,x)  dT  =  number  of  PKA's  of  energy  dT  about  T  produced 
per  unit  time  per  unit  bombarding  beam  flux 
as  an  ion  beam  with  distribution  0  (E,X)  traverses 
a  perpendicular  unit  area  located  at  depth  x.  (4) 

We  have 

W(T,x)  =  /dE  0(E,x)  K12(E,T)  dR12(E)/dRp.  (5) 


Energy  Deposition  by  a  PKA 

Consider  a  PKA  created  at  the  point  A  at  depth  Xj  with  energy  Tj  and 
velocity  direction  v..  We  will  assume  no  range  straggling  for  Z2  atoms.  Our 
model  has  range  and  energy  straggling  only  for  the  beam  (Z^)  atoms.  We  also 
assume  that  the  PKA  travels  in  a  straight  line,  coming  to  rest  at  a  point  C. 
Let  B  be  an  intermediate  point  at  depth  x2,  where  the  PKA  has  energy  T2> 
Let  n  be  the  cosine  of  the  angle  between  Vj  and  the  x-axis: 

H  (Ej,Tj)  =  Vj  .  x  ;  (6) 


the  kinematics  of  elastic  collision  determines  /j  as  a  function  of  Ej  and  Tj, 
where  E{  is  the  energy  of  the  Zj  beam  ion  creating  the  PKA.  We  have 

x2-Xj+  n  (Ej,Tj)  Rp22  (T2)  =  Rp22  ^) 

This  equation  determines  T2  as  a  function  of  T j  and  Ej: 

T  -  T  (E  T  )  (8) 

2  _  ‘2  vtr  V 


Now  let  Sp(x2,  x ,  ,T , )  dx2  be  the  damage  energy  deposited  by  this  PKA  as  it 
traverses  a  aepth  ax2  around  the  point  B.  We  have 

Sp(x2,Xj,  Tj)  =  Sp(x2-x1,T][) 

=  /l/d(E1,T1)^i722(T2)Sm22(T2)  ^dR22(T2)/dRpj,  (9) 

where  r]  9(T)  is  the  energy  partition  factor  of  Lindhard,  Nielsen,  Scharff,  and 
Thompson  11). 
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Energy  Deposition  by  the  Bombarding  Beam 


For  a  bombarding  beam,  the  damage  energy  deposited  at  the  point  x 2 
is  obtained  by  multiplying  the  expression  (9)  by  the  number  of  PKA's  of  energy 
T.  produced  at  x^  and  summing  over  Tj  and  x^.  We  will  denote  by 
tne  damage  energy  including  recoil  effects: 


SDR(x2)  = 


^2  x 

mhn 


Sp(x2,Xi ,  T j)W  (Tj ,  X|). 


Equation  (5)  gives 


X.  T 

2  in 


SDR  =  /  d*i  /  dTi  j dEj  Sp  (x2,  Xj,  Tj^Ej,  Xj) 
o  E 

K12(E1’T1)  ^Rj^EjJ/dRp^ 


Let  Xj  be  the  stopping  depth  of  the  Zj  beam  ion  which  creates  the  PKA  of 
energy  T  at  depth  x ^ ;  this  Zj  ion  has  energy  E^  at  the  depth  Xj.  Using  the 
methods  of  Reference  (1 1),  we  can  employ  the  range  distribution  f(x^)  of  Zj 
ions  to  estimate  the  flux  0  (E. ,  x.)  of  equation  (1 1),  obtaining 

X„  T  1 

2  m  oo 


SDR  (x2*  " 


/<h>/dT'/ 


dx^  Sp  (x2-Xj,  T|) 


K12  (EpTj)  (dR12/dRp)  f(x3) 


Now,  interchange  the  order  of  integration  of  x^  and  Tj,  and  define 


SM(x2'xl,El)  = 


=  /dT‘K 


12  ^1’  ^p  Cx2_x i f) 


with  Sp  given  by  equation  (9). 


S|^dx2  is  the  damage  energy  dE^  deposited  at  dx2  around  x 2  by  all 
the  PKA's  created  at  depth  Xj  by  a  single  beam  ion  Zj  with  energy  Ep 


Define  SM2  (x^Xj)  as 


S  (x 
^MS  vx2 


-Xj)  =  f  dx3  SM(x2-Xl,  EjJ^dRj^E^/dRp^ffxj) 
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!>MS  is  the  damage  energy  function  for  a  penetrating  beam  with  the  energy 
straggling  given  by  the  distribution  function  0(Ej,Xj).  We  now  get  the  desired 
energy  profile  S^R  of  equation  (12)  as 


SDR(x2) 


■/ 


dxl  SMS(x2"xl)* 
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II.  HEAVY  ION  AND  NEUTRON  DAMAGE  STUDIES 


A.  The  Effect  of  Alloy  Additions  on  the  Swelling  of  Neutron- 
Irradiated  Nickel 


(3.  A.  Sprague,  K.  B.  Roarty,  and  F.  A.  Smidt,  Jr. 
Thermostructural  Materials  Branch, 

Material  Science  <3c  Technology  Division) 


Background 


Swelling  (or  volume  change)  of  materials  caused  by  high-fluence  neutron 
irradiation  at  elevated  temperatures  is  a  major  engineering  design  problem 
that  limits  the  performance  of  advanced  nuclear  reactors.  This  swelling  is 
caused  by  the  agglomeration  of  radiation-produced  vacancies  into  microscopic 
cavities  or  voids.  The  design  of  new  low-swelling  alloys  for  nuclear  applica¬ 
tions  requires  fundamental  knowledge  of  the  effect  of  individual  alloying 
elements  in  simple  binary  alloy  systems,  since  individual  effects  cannot  be 
readily  studied  in  the  more  complex  commercial  alloys.  The  present  study  was 
undertaken  to  examine  the  effects  of  small  additions  of  aluminum,  molyb¬ 
denum,  silicon,  and  titanium  on  void  formation  in  neutron-irradiated  nickel. 
These  data  were  then  compared  with  the  results  of  previous  work  (12)  that 
involved  nickel-ion  irradiation  of  the  same  set  of  alloys.  This  comparison 
provides  insight  into  the  effects  of  alloying  on  both  the  swelling  of  nickel  and 
the  correlation  between  neutron  and  heavy-ion  damage. 


Progress 


The  pure  nickel  control  material  and  the  nickel  alloys  used  in  the  present 
study  were  taken  from  the  alloy  series  previously  examined  following  nickel- 
ion  irradiation  (12).  The  alloys  Ni-L  at.%Al,  Ni-1  at.%Mo,  Ni-1  at.%Si,  and  Ni- 
1  at.%Ti  were  selected  for  the  first  examinations  following  neutron  irradia¬ 
tion.  The^ion  irradiation  experiments  had  shown  that  the  swelling  in  Ni  (at  8.1 
dpa*  625  C)  was  reduced  by  a  factor  of  40  by  the  A1  addition,  and  that 
significant  void  formation  was  completely  suppressed  by  the  Mo,  Si,  and  Ti 
additions. 


The  alloys  and  the  control  material  were  all  prepared  by  arc-melting 
''■'65-gram  buttons  in  an  argon  atmosphere,  which  had  been  purified  by 
titanium  gettering.  After  several  remelts,  a  3-mm  thick  slice  was  sectioned 
from  the  center  of  each  button.  The  slices  were  cold-rolled  to  0.1 -mm  thick 
foil,  from  which  3-mm  diameter  disks  were  punched.  With  the  exception  of  the 
A1  alloy,  these  specimens  were  annealed  for  one  hour  at  990°C  in  a  1  x  10”6 
torr  vacuum.  The  Ni-1%A1  specimens  were  similarly  annealed,  but  the 
temperature  was  reduced  to  800°C  to  reduce  the  loss  of  A1  from  the  foils. 


The  original  reference  calculated  dpa  values  using  a  displacement  energy 
of  25  eV  for  nickel,  while  the  present  report  uses  40  eV  to  conform  with 
current  convention. 
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I 


The  disks  were  irradiated  in  the  Experimental  Breeder  Reactor  No.  II  (EBR-II)  fast 
reactor  in  Idaho  Falls,  Idaho,  to  a  neutron  fluence  of  2.2  x  10 22  neutrons/cm  , 
E>  0.1  MeV,  at  a  temperature  of  455°C.  Following  irradiation,  disks  were 
thinned  by  dual-jet  electropolishing  in  a  solution  of  250  ml  methyl  alcohol,  150 
ml  n-butyl  alcohol,  and  10  ml  perchloric  acid  at  a  temperature  of  -60°C  and  a 
voltage  of  35  V.  Void  and  dislocation  structures  were  studied  by  transmission 
electron  microscopy. 

The  dislocation  structures  in  all  of  the  materials  consisted  primarily  of 
tangles  of  line  dislocations.  This  is  illustrated  in  Fig.  1,  which  shows  the 
structure  observed  in  the  pure  Ni.  As  is  normally  observed  in  specimens 
containing  voids,  the  dislocation  lines  tended  to  interconnect  the  voids. 

The  qualitative  features  of  the  void  distributions  observed  in  Ni,  Ni- 
1%A1,  Ni-l%Mo,  Ni-l%Si,  and  Ni-1  %Ti  are  shown  in  Figs.  2-6,  respectively. 
Quantitative  parameters  measured  from  these  void  structures  are  given  in 
Table  1.  The  following  features  of  void  formation  in  these  materials  can  be 
seen: 

1.  The  void  size  distribution  in  the  pure  Ni  was  very  broad,  and  the 
voids  were  formed  heterogeneously,  with  clusters  of  small  voids 
dispersed  among  larger  ones. 

2.  In  the  Ni-Al  alloy,  the  void  size  distribution  was  narrower,  and  the 
spatial  distribution  of  the  voids  was  more  homogeneous,  but  the 
average  size,  number  density,  and  swelling  were  similar  to  the  pure 
Ni. 

3.  The  Mo  and  Ti  additions  both  increased  the  swelling  by  approxi¬ 
mately  a  factor  of  two,  but  by  different  mechanisms.  The  Mo 
increased  the  void  nucleation  but  did  not  affect  the  void  growth, 
leading  to  a -higher  density  of  similarly-sized  voids.  The  Ti,  on  the 
other  hand,  did  not  appear  to  alter  the  void  nucleation,  but 
increased  void  growth,  leading  to  a  similar  density  of  larger  voids. 

4.  The  Si  addition  produced  the  only  reduction  in  swelling  of  the  four 
alloying  elements,  decreasing  it  by  approximately  a  factor  of  ten. 
The  Si  did  not  significantly  change  the  void  number  density,  but 
considerably  reduced  the  mean  void  size. 

These  observed  alloying  effects  on  void  formation  differ  significantly 
from  the  previously  reported  nickel-ion  e  -'eriments  on  the  same  alloys  (12). 
As  noted  above,  all  four  of  the  alloying  fc.oments  reduced  or  completely  sup¬ 
pressed  the  swelling  in  the  Ni-ion-irradiated  Ni,  by  reducing  the  void  nuclea¬ 
tion  (assuming  that  the  ion-irradiated  specimens  did  not  contain  a  high  density 
of  voids  smaller  than  the  ^1.5  nm  resolution  limit  of  the  microscope).  A 
thorough  understanding  of  the  differences  between  neutron  and  nickel-ion 
irradiations  of  these  alloys  would  require  a  much  more  extensive  investigation 
of  irradiation  temperatures  and  damage  levels  for  both  irradiating  particles. 
Experimental  and  theoretical  models  which  are  currently  available,  however, 
offer  some  potential  explanations.  A  study  of  the  swelling  produced  in  piye  Ni 


The  dislocation  structure  observed_in  pure  nickel  following  irrad 
tion  with  2.2  x  10  neutrons/cm  ,  E>0.1  MeV,  at  455°C.  1 

dislocation  structures  observed  in  Ni-1  at.%Al,  Ni-1  at.%Mo,  Ni 
at.%Si,  and  Ni-1  at.%Ti  were  qualitatively  similar. 


Fig-  2 


Vo^  distribution  observed  in  pure  Ni  following  irradiation  with  2.2  x 
10  neutrons/cm  ,  E  >0. 1  MeV,  at  *f55°C. 


15 


Fig.  3 

Void  distribution  observed  in  Ni-1  at.%Al  following  iradiation  with 
2.2  x  10  neutrons/cm  ,  E  >0.1  MeV,  at  4 53°C. 
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Void  distribution  observed  in  Ni-I  at.%Mo  following  irradiation  with 
2.2  x  10Z/  neutrons/cm  ,  E  >0.1  MeV,  at  455°C. 
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Fig.  5 

Void  distribution  observed  in  Ni-1  at.%Si  following  irradiation  with 
2.2  x  lCr^  neutrons/cm  ,  E  >0.1  MeV,  at  455°C. 


Fig.  6 

Void  distribution  observed  in  Ni-1  at.%Ti  following  irradiation  with 
2.2  x  1CT  neutrons/cm  ,  E  >0.1  MeV,  at  455°C. 
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table  1 

VOID  FORMATION  IN  NICKEL  ALLOYS 
(2.2  x  1022  n/cm2,  455°C) 


Material 

Void  Density 
(cm'7) 

Mean  Void 

Dia  (nm) 

Swelling 

(%) 

Ni  Control 

Ni-1%  A1 

Ni-1%  Mo 

Ni-1%  Si 

Ni-1%  Ti 

3.5  x  lo\l 

4.2  x  10  ,, 

1.33  xl^5 

3.8  x  10}7 

2.8  x  1014 

43.5 

48.8 

41.0 

21.0 

62.1 

2.56 

3.04 

5.64 

0.27 

4.91 

by  neutron  and  Ni-ion  irradiations  (13)  concluded  that  different  void  nucleation 
mechanisms  were  controlling  void  formation  in  the  two  cases.  In  the  ion 
irradiation  experiments,  the  displacement  rate  was  4.4  x  10~z  dpa/sec,#  which 
produces  large  super  saturations  of  vacancies  and  interstitials,  which  can  lead 
to  significant  void  nucleation  in  the  absence  of  insoluble  gag>es  (14).  In  a  fast 
reactor  irradiation,  however,  displacement  rate  is  M  x  10"°  dpa/sec,  and  the 
helium  produced  by  transmutation  reactions  is  probably  responsible  for  driving 
the  void  nucleation.  It  is  therefore  possible  that  all  four  alloying  elements  in 
this  study  could  greatly  reduce  the  homogeneous  nucleation  rate  in  a  Ni-ion 
irradiation  while  either  having  no  effect  on  the  helium-assisted  nucleation 
mechanism  that  is  dominant  in  a  neutron  irradiation,  or  even  enhancing  this 
mechanism.  A  second  factor  which  could  have  influenced  these  experiments  is 
radiation-induced  segregation  of  the  alloying  elements.  In  other  studies  (15), 
it  has  been  found  that  alloying  elements  which  trap  vacancies  or  interstitials 
will  not  only  affect  the  void  nucleation  and  growth,  but  will  also  segregate 
with  respect  to  internal  point  defect  sinks,  such  as  voids  and  dislocations.  This 
segregation  can  change  the  rates  at  which  these  sinks  absorb  vacancies  and 
interstitials,  as  well  as  reduce  the  amount  of  a  solute  element  in  solution  in  an 
alloy  matrix.  Both  of  these  effects  can  alter  the  void  formation  in  an  alloy. 
The  extent  of  radiation-induced  solute  segegation  is  a  strong  function  of  the 
damage  rate,  so  that  it  ^vill  be  very  different  in  neutron  and  nickel-ion 
irradiations. 

The  present  neutron-irradiation  experiments,  together  with  the  previous 
ion-irradiation  work  on  the  same  alloys,  have  demonstrated  that  low  concen¬ 
trations  of  solid-solution  alloying  elements  can  significantly  change  the  void 
formation  behavior  of  nickel.  More  importantly,  these  studies  indicate  that 
caution  should  be  exercised  in  evaluating  high-dose-rate  ion  irradiations  of 
complex  alloys.  The  effects  of  damage  rate  on  all  important  processes, 
including  nucleation  mechanisms  and  solute  segregation,  must  be  considered  if 
such  ion  irradiation  experiments  are  to  accurately  reflect  the  relative  swelling 
behavior  of  the  alloys  during  reactor  irradiations. 
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B.  The  Microstructure  of  Neutron  Irradiated  Refractory 

Metals  and  Alloys 

(3.  A.  Sprague,  F.  A.  Smidt,  3r.,  and  3.  R.  Reed, 

Thermostructural  Materials  Branch, 

Material  Science  &  Technology  Division) 

Background 

Alloys  based  on  the  bcc  refractory  metals  are  of  interest  for  high 
temperature  nuclear  systems,  such  as  fusion  power  reactors,  because  of  the  high 
operating  temperatures,  and  thus  high  power  densities,  possible  with  this  class  of 
materials.  The  engineering  application  of  refractory  alloys  will  require  the 
expansion  of  our  current  data  base  on  these  materials  and  their  behavior  during 
high-fluence  irradiation  at  elevated  temperature.  The  present  experiments,  which 
are  closely  related  to  the  work  of  the  CORES  program,  were  performed  as  part  of 
the  Fusion  Materials  Program  of  the  Office  of  Fusion  Energy,  U.  S.  Department  of 
Energy.  The  objective  was  to  study  the  microstructural  response  of  molybdenum, 
TZM,  niobium,  niobium-1%  zirconium,  vanadium,  and  vanadium-20%  titanium  to 
high  fluence  neutron-irradiation  at  650°C.  These  results  supplement  previous  work 
by  a  number  of  investigators  to  provide  a  basis  for  <he  more  systematic  and 
comprehensive  experiments  that  will  be  required  to  adequately  define  the 
irradiation  response  of  refractory  metals  and  alloys. 

Progress 

Experimental  Procedures.  The  pure  Mo  (>99.95%),  Nb  (>99.99%),  and  V 
(>  99.95%)  were  obtained  from  Materials  Research  Corporation  in  the  form  of 
50-/im  thick  foil.  (Note  that  the  above  nominal  purities  are  for  substitutional 
solutes  only.)  The  TZM  (Mo-0.5  wt.%  Ti-0.1  wt.%  Zr),  Nb-1  wt.%  Zr,  and  V-20  wt.% 
Ti  materials  were  rolled  to  50-fim  thick  foil  from  commercial-grade  stock. 
Following  the  punching  of  3-mm  dia.  disks,  the  metals  and  alloys  were  all  annealed 
in  a  vacuum  of  1.3  x  10"^  Pa  for  1  hour  at  1430°C,  except  for  the  V  and  V-20  Ti, 
which  were  annealed  at  885°C.  The  foils  were  packed  into  close-fitting  stainless 
steel  tubes  which  were  evacuated,  back-filled  with  helium,  and  welded  shut.  The 
foil  tubes  were  loaded  into  sodium-filled  heat  pipe-controlled  irradiation  capsules 
(16)  for  irradiation  in  the  EBR-II  reactor.  Two  capsules  were  used  in  the  present 
experiments,  both  designed  to  operate  at  650  +10°C.  The  first  (EBR-II  subassembly 
X-200),  which  contained_all  materials  except  the  V-20  Ti,  was  irradiated  in  Row  7 
to  a  fluence  of  5.4  x  KrZ  neutrons/cm  ,  E>0.1  MeV.  The  second  capsule  (EBR-II 
subassembly  X-255)  containing  the  Mo  and  V-20  _^i  was  irradiated  in  several 
positions,  mostly  in  Row  2,  to  a  fluence  of  8.4  x  10^Z  neutrons/cm  ,  E>0.1  MeV. 
Assuming  a  40  eV  displacement  energy  for  all  the  materials,  the  lower  neutron 
fluence  corresponds  to  22  dpa  for  Mo  and  TZM,  21  dpa  for  Nb  and  Nb-1  Zr,  and  24 
dpa  for  V.  The  higher  neutron  fluence  corresponds  to  35  dpa  for  Mo  and  38  dpa  for 
V-20Ti.  Following  irradiation,  the  foils  were  electrochemically  thinned  and 
examined  by  transmission  electron  microscopy. 
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Results  and  Discussion 


The  experimental  results  will  be  discussed  in  each  of  the  following 
subsections  by  comparing  each  pair  of  pure  metal  and  alloy,  in  turn.  The 
quantitative  void  distribution  data  for  those  materials  that  exhibited  significant 
void  formation  are  all  collected  in  Table  2,  and  will  be  referenced  in  each 
subsection,  as  needed. 


Mo  and  TZM.  The  microstructures  of  the  Mo  at  both  fluences  and  the  TZM  at  the 
lower  fluence  were  similar,  consisting  of  imperfectly  ordered  void  lattices  and 
moderate  densities  of  dislocation  loops  and  line  segments.  Examples  of  these 
microstructures  are  shown  in  Fig.  7.  The  total  dislocation  density  in  each  of  the 
specimens  examined  was  in  the  range  5  -  7  x  1CT  cm”  .  A  low  density  of  globular 
precipitates  was  observed  in  the  TZM.  In  the  pure  Mo,  a  void  lattice  parameter  of 
30  nm  was  measured  for  both  fluences,  which  is  in  good  agreement  with  previous 
lower-fluence  measurements  of  void  lattice  parameter  as  a  function  of  irradiation 
temperature  (17,18).  Some  ordering,  on  approximately  the  same  scale  as  in  the  Mo, 
could  be  detected  in  the  TZM,  but  a  good  void  lattice  parameter  determination 
could  not  be  made,  due  to  the  relatively  large  size  of  the  voids  and  the 
imperfections  in  the  ordering.  Referring  to  the  void  distribution  data  in  Table  2,  it 
can  be  seen  that  the  swelling  measured  in  TZM  at  the  lower  fluence  was  slightly 
greater  than  that  in  the  pure  Mo,  although  the  difference  was  within  the  error 
limits  for  these  determinations.  At  the  higher  fluence,  where  the  only  data  were 
for  the  pure  Mo,  the  swelling  was  actually  smaller  than  at  the  lower  fluence, 
although  the  difference  was  again  within  error  limits.  Both  void  densities  and  mean 
diameters  in  the  Mo  were  similar  at  both  fluences.  The  TZM  contained  a  lower 
density  of  significantly  larger  voids,  a  number  of  which  appeared  to  have  coalesced 
from  multiple  smaller  voids  formed  earlier  in  the  irradiation.  A  final  observation 
was  that  the  voids  in  the  pure  Mo  were  surrounded  by  the  strong  strain  fields, 
which  were  easily  visible  in  dark  field,  using  a  <  110>g-vector.  No  strain  fields  of 
similar  magnitude  were  observed  surrounding  the  voids  in  TZM. 


The  present  observation  of  the  relative  swelling  in  Mo  and  TZM  agrees 
fairly  well  with  previously  reported  results  of  other  investigators,  although  there 
are  significant  differences  in  relative  void  sizes  and  number  densities.  At  low 
neutron  fluences  (M0ZU  n/cml,  Bentley,  et.  al  (19)  reported  significant  void 
formation  in  Mo  over  the  entire  temperature  range  330  -  850°C,  while  they 
observed  significant  void  formation  in  TZM  only_£pr  "intermediate"  temperatures  in 
this  range.  At  higher  fluences  of  2.5  -  4.4  x  10  z  n/cm  ,  E>0.1  MeV,  Bentley  and 
Wiffen  (17)  observed  comparable  swelling  in  Mo  and  TZM  at  425  and  585  C,  and 
slightly  greater  swelling  in  TZM  at  790  and  1000°C.  The  void  number  densities 
observed  in  pure  Mo  in  the  present  study  are  consistent  with  the  temperature 
dependence  of  void  number  density  reported  by  Bentley  and  Wiffen.  There  is  a 
significant  difference,  however,  between  the  void  population  reported  by  Bentley 
and  Wiffen  in  TZM  and  those  observed  in  the  present  study.  Compared  with  pure 
Mo,  they  observed  a  higher  density  of  smaller  or  equal  size  voids  in  TZM  at  all 
temperatures  between  425  and  1000°C,  while  the  present  study  found  a  lower 
density  of  larger  voids  in  TZM.  Another  difference  between  the  present  results  and 
those  in  the  literature  is  with-Jhose  reported  by  Sikka  and  Moteff  (18)  for 
irradiation  to  a  fluence  of  I  x  n/cnrr  Ei  I  MeV,  at  temperatures  between  430 
and  1000°C.  Although  this  fluence  corresponds  to  21  dpa,  the  same  as  the  present 
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TABLE  2 


VOID  DISTRIBUTIONS  OBSERVED  AFTER  NEUTRON  IRRADIATION  AT  650°C 


Material 


Fluen^e 
(n/cm  , 
E>  0.1  MeV) 


Void  Density 
(cm'3) 


Mean  Diam. 
(nm) 


Swelling 

(%) 


Representative  void  structures  in  irradiated  Mo  and  T2M:  (a)  Mo, 
imaged  with  a  [ill]  beam  direction  showing  void  ordering;  the  voids 
appear  dark  due  to  a  combination  of  overfocus  and  strain  contrasts; 
(b)  TZM,  in  which  several  coalesced  voids  can  be  seen. 


study,  their  reported  temperature  dependence  implies  only  approximately  10 
percent  of  the  presently  observed  swelling.  The  reasons  for  the  above  discre¬ 
pancies  are  not  known  at  this  time. 

Nb  and  Nb  -  1%  Zr.  The  Nb  contained  a  very  high  density  of  small  voids,  as 
illustrated  in  Fig.  8a,  as  well  as  a  much  lower  density  of  small  ('MO  nm) 
precipitates  in  the  matrix.  The  void  lattice  parameter,  assuming  a  bcc  structure, 
was  24  nm.  The  swelling,  as  seen  from  the  data  in  Table  2,  was  considerably  less 
than  that  in  the  other  pure  metals  or  the  TZM.  Prior  to  irradiation,  the  Nb-1  Zr 
contained  a  few  precipitates  at  the  grain  boundaries.  After  irradiation,  the  alloy 
contained  additional  precipitates,  both  at  the  grain  boundaries  and  in  the  matrix,  as 
seen  in  Fig.  8b.  A  few  voids,  ^5  nm  in  diameter,  were  observed  in  isolated  areas  of 
the  foil,  but  these  did  not  represent  any  significant  swelling.  This  observation  of 
suppression  of  swelling  in  Nb-lZr  agrees  with  the  results  of  Michel  and  Moteff  (20) 
and  Tang  and  Moteff  (21),  who  both  reported  significant  swelling  in  Nb-lZr  ordx 
over  ^he  temperatur^j-ange  700  -  900°C  at  a  neutron  fluences  of  2.5  x  KT 
n/cm  (20)  and  5  x  10  n/cnn  (21),  E  >0.1  MeV.  The  void  distribution  observed  in 
Nb  in  the  present  study,  however,  was  quite  different  from  that  reported  by  Michel 
and  Moteff  in  that  the  voids  observed  here  were  smaller,  and  the  number  density 
was  approximately  an  order  of  magnitude  larger.  The  reasons  for  this  difference 
cannot  be  determined  unambiguously,  but  they  are  probably  related  to  the  impurity 
contents  of  the  materials  in  the  two  studies.  First,  the  material  used  in  Michel  and 
Moteff's  work  was  commercial-grade  Nb,  as  opposed  to  the  high-purity  stock  used 
in  the  present  experiment.  Second,  the  specimens  in  the  present  study  were 
annealed  as  50-Um  thick  foil  in  a  1.3  x  1(j  Pa  vacuum  at  1430°C.  Although  the 
specimen  holder  in  the  furnace  was  made  of  tantalum,  considerable  contamination 
(on  the  order  of  a  few  thousand  parts  per  million)  by  interstitial  impurities  probably 
occurred  (22).  This  conclusion  is  reinforced  by  the  observed  precipitation.  The 
preirradiation  annealing  conditions  for  the  material  used  in  Michel  and  Moteff's 
work  were  not  specified.  It  seems  reasonable,  therefore,  that  differences  in 
substitutional  and  interstitial  impurity  contents  were  responsible  for  the  different 
void  distributions. 

V  and  V-20%  Ti.  The  microstructure  of  the  irradiated  pure  V  was  quite  different 
from  those  of  the  Mo  and  Nb,  as  shown  in  Fig.  9a.  The  preirradiation  annealing 
treatment  had  produced  a  fine  grain  size  (5  -  10 pm),  and  after  irradiation,  each 
grain  contained  a  moderate  density  of  fairly  large  voids.  The  voids  were  denuded 
near  eaqji  graip  boundary.  The  dislocation  density  in  the  grain  interiors  was  low, 
^  1  x  10*  cm  ,  although  many  dislocations  were  observed  in  subgrain  boundaries. 
No  precipitates  were  observed.  The  preirradiation  microstructure  of  the  V-20Ti 
contained  a  heterogeneous  distribution  of  precipitates,  the  diffraction  patterns  of 
which  were  consistent  with  hexagonal  Ti02  (rutile).  The  maximum  observed  local 
number  density  of  these  precipitates  was  2  xiD  J  cm"  ,  while  the  average  density 
over  a  larger  specimen  volume  was  ^4  x  101*  cm"  .  A  detailed  characterization 
of  the  precipitate  size  distributions  was  not  made,  but  two  morphologies  were 
observed:  blocky  precipitates  with  a  150  nm  mean  diameter,  and  rod-shaped,  MOO 
nm  in  crgp  section  and ^00  -  400  nm  long.  Irradiation  of  this  material  at  650°C  to 
8.4  x  10 neutrons/cm  ,  E  >0.1  MeV,  produced  no  observable  voids,  but  did  result 
in  significant  growth  of  the  precipitates.  The  mean  diameter  of  the  blocky 
particles  increased  to  500  nm;  the  mean  cross  section  of  the  rod-shaped  precipi¬ 
tates  was  150  nm,  while  the  maximum  length  of  the  rods  increased  to  1500  nm.  A 
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Nb,  containing  a  high  density  of  small  voids;  (b)  Nb-1%  Zr,  with  no 
void  formation,  but  extensive  matrix  and  grain  boundary  precipi¬ 
tation.  Note  the  different  magnifications  in  the  two  micrographs. 


Kepresentative  microstructures  in  irradiated  V  and  V-20%  Ti:  (a)  V, 
showing  the  denuding  of  voids  at  boundaries  of  the  small  grains;  (b) 
V  -  20%  Ti,  containing  elongated  TiC^  precipitates  which  increased 
in  maximum  length  from  400  to  1400  nm  during  irradiation. 


cluster  of  the  elongated  precipitates  in  the  V-20Ti  is  shown  in  Fig.  9b. 


The  complete  suppression  of  void  formation  by  the  addition  of  20%  Ti  to  V 
agrees  with  the  results  of  B^itley  apd  Wiffen  (23),  who  examined  V-20Ti  irradiated 
to  f^uences  of  1.3  -  6.1  x  10z  ^  n/cm  ,  E  >0.1  MeV,  at  temperatures  between  470  and 
780°C.  Before  irradiation,  their  specimens  contained  precipitates  similar  in  size, 
shape,  and  maximum  number  density  to  those  observed  prior  to  irradiation  in  the 
present  study,  although  they  reported  that  the  specimen  regions  containing  this 
maximum  number  density  represented  only'v^  of  the  specimen  volume.  They 
observed  no  change  in  this  preirradiation  precipitate  structure  after  irradiation. 
Carlander  (24)  reported  a  few  voids  in  V-20Ti  after  irradiation  to  6  x  l(j  n/cm  , 
E  >0.1  MeV,  at  600°C,  although  these  voids  did  not  represent  any  significant 
swelling.  He  did  not  report  on  the  precipitation  behavior  of  this  alloy,  but  a 
number  of  apparent  rod-shaped  precipitates  can  be  seen  in  his  published  micro¬ 
graph  (Fig.  3  of  Carlander,  Reference  24). 


The  void  microstructures  in  the  pure  V  are  difficult  to  correlate  with  those 
reported  previously  for  lower  fluence  irradiations  (24,25),  due  to  differences  in 
starting  materials  2^nd  irradiation  temperatures,  but  the  observation  of  3.6% 
swelling  at  5.4  x  1(j  n/cnn  is  consistent  with  approximately  linear  swelling  with 
fluence,  when  compared- with  the  1.54%  swelling  reported  by  Carlander  (24)  for 
irradiation  to  2.5  x  10Z  n/cnn  at  6Zi°C  and  the  1.47%  swelling  reported  by 
Stiegler  (25)  for  irradiation  to  1.4  x  10  Z  n/cnn  at  600°C  (all  fluences  for  E  >0.1 
MeV). 


Conclusion 


The  present  microstructural  observations  of  refractory  metals  and  alloys 
after  high-fluence  neutron  irradiation  at  650°C  indicate  that,  at  least  at  this 
temperature,  void  swelling  will  probably  not  be  the  most  serious  obstacle  to  the 
application  of  this  class  of  alloy.  The  swelling  of  even  the  pure  metals  is  only 
moderate,  and  with  the  possible  exception  of  Mo  alloys,  it  appears  possible  to 
control  void  formation  further  by  alloying.  Changes  in  mechanical  properties  due  to 
precipitation  of  interstitial  impurities  and  the  formation  of  high  densities  of  small 
voids  may  be  far  more  serious  problems.  Theoretically,  the  precipitation  could  be 
controlled  by  reducing  the  concentration  of  interstitial  contaminants,  but  more 
comprehensive  studies  on  thoroughly  characterized  materials  will  be  required  to 
examine  this  possibility.  If  detrimental  changes  in  microstructure  can  be  avoided 
by  careful  control  of  impurities,  fabrication  technologies  and  reactor  operating 
conditions  will  still  have  to  be  developed  to  maintain  this  impurity  control  in  large 
structures.  If  these  developments  can  be  made,  however,  refractory  metal  alloys 
appear  to  be  potentially  valuable  for  high  temperature  nuclear  applications. 
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III. 


LIGHT  ION  BEAM  STUDIES 


A.  Deuterium  Retention  in  Aluminum  and  Stainless  Steel  With  and 
Without  Helium  Loading 

(L.  A.  Beach,  Cyclotron  Applications  Branch, 

Radiation  Technology  Division) 

Fast  neutrons  in  high  flux  reactors  generate  hydrogen  and  helium  in  reactor 
structural  materials  by  (n,p)  and  (n,  a)  reactions.  Helium  is  known  to  accumulate 
in  metals  and  alloys  under  reactor  operating  conditions  and  its  accumulation  is 
intimately  related  to  swelling  and  other  changes  in  mechanical  properties.  It 
has  generally  been  assumed  that  hydrogen,  although  generated  in  larger  amounts 
than  helium,  diffuses  or  permeates  through  the  material  and  consequently  not 
causing  any  mechanical  property  changes.  However,  retention  of  appreciable 
amounts  of  hydrogen  could  lead  to  deleterious  effects. 

Lanford  (26)  reported  that  hydrogen  implanted  into  single  crystal  silicon 
wafers  was  retained  indefinitely  but  that  hydrogen  retention  in  tungsten  was 
enhanced  with  helium  predamage  with  little  retention  in  samples  of  tungsten 
without  predamage.  Nowicki  (27),  based  upon  hydrogen  profile  measurements  by 
Pieper  at  NRL,  showed  that  large  amounts  of  hydrogen  were  present  in  rf- 
spattered  A^O^  films  deposited  by  planer  magnetron,  it  remained  stable  and 
affects  mechanical  properties.  Picraux  (28)  has  shown  enhanced  hydrogen  and 
deuterium  retention  in  molybdenum  with  helium  predamage,  both  for  single 
crystal  samples  and  rolled  foil.  Since  helium  is  known  to  become  substitutional 
by  falling  into  vacancies  created  by  atomic  displacements  and  positrons  have 
been  observed  to  attach  to  these  helium  filled  vacancies,  the  possibility  exists 
that  hydrogen,  which  would  normally  remain  interstitial  and  diffuse  through  the 
metal,  may  attach  itself  to  these  filled  vacancies. 

In  controlled  thermonuclear  systems  high  fluxes  of  14-MeV  neutrons  will 
bombard  structural  materials  including  first  walls,  blankets  and  cyrogenic 
cooling  systems.  The  generation  of  hydrogen  and  helium  will  be  enhanced  by  the 
higher  cross-sections  and  in  addition  large  fluxes  of  deuterium  and  tritium 
escaping  from  the  plasma  will  strike  the  first  walls.  The  question  of  hydrogen 
retention  being  enhanced  by  helium  loading  or  radiation  damage  is  of  impor¬ 
tance,  particularly  at  higher  temperatures.  Of  course,  the  payoff  will  be  if 
retained  hydrogen  causes  changes  in  mechanical  properties. 

As  a  first  step,  samples  of  aluminum,  stainless  steel,  molybdenum  and 
tantalum  were  implanted  with  deuterium  both  with  and  without  prior  helium 
implantations.  For  the  helium  loading,  a  60  MeV  a  -particle  beam  from  the 
cyclotron  was  degraded  sufficiently  to  stop  within  the  sample  and  due  to 
straggling  give  a  broad  profile.  A  helium  implant  in  a  1.2  mil  aluminum  foil  was 
degraded  to  give  a  profile  peaking  at  a  depth  of  0.5  mils  as  shown  by  the  solid 
curve  in  Fig.  10.  The  degradation  of  the  alpha  beam  was  selected  to  give  a 
profile  peaking  at  0.17  mils  in  a  0.6  mil  stainless  steel  sample,  see  Fig.  11.  A 
proton  beam  was  not  used  as  the  hydrogen  isotope  implant  because  previous 
measurements  (29)  showed  surface  contamination  of  the  order  of  micrograms  of 
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Fig.  10 

Helium  profile  in  a  1.2  mil  aluminum  sample  implanted  with  helium 
and  deuterium.  The  solid  curve  is  the  predicted  implantation 
obtained  from  the  computer  code  normalized  by  integrated  cur¬ 
rent.  The  data  points  are  from  a  -  a  scattering  when  a  36.5  MeV 
alpha  beam  impinges  on  the  sample. 
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Fig.  11 

Helium  profile  in  a  0.6  mil  stainless  steel  sample  implanted  with 
both  helium  and  deuterium  as  in  Fig.  1.  The  solid  curve  is  the 
predicted  helium  profile. 
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hydrogen  per  cm  complicate  hydrogen  profiling  by  p-p.  elastic  scattering.  A  22.5 
MeV  deuteron  beam  was  similarly  degraded  to  give  deuteron  profiles  peaking  at 
the  same  depths  in  aluminum  and  stainless  steel.  The  molybdenum  (2.0  mils)  and 
tantalum  (1.0  mil)  were  similarly  loaded  but  these  specimens  were  too  thick  to 
profile  at  the  loadings.  The  predicted  profiles  are  obtained  from  a  computer  code 
written  by  Pieper  that  gives  implantation  concentrations  in  atomic  parts  per 
million  (appm)  per  microcouloumb  of  implantation  beam. 

The  loadings  of  two  aluminum  and  two  stainless  steel  specimens  are  shown 
in  Table  3.  Helium  and  deuterium  profiles  were  measured  by  scattering  with  a 
36.5  MeV  alpha  beam  and  by  d-d  scattering  of  a  22.5  MeV  deuteron  beam, 
respectively.  A  measured  helium  profile  of  the  A1  sample  implanted  with  both 
a's  and  d's  is  shown  in  Fig.  10  in  reasonable  agreement  with  the  predicted  profile 
normalized  by  the  integrated  current  during  helium  implantation.  Due  to  the  low 
implantation  levels,  statistics  are  poor  and  accidental  backgrounds  from  host 
scatteirngs  generated  problems.  The  deuteron  profile  for  this  same  aluminum 
sample  is  shown  in  Fig.  12  along  with  the  predicted  profile.  The  deuteron  profile 
for  the  aluminum  sample  not  preloaded  with  helium  is  similar  and  agrees  well 
with  the  predicted  loading  indicating  that  the  deuterium  is  being  retained  in 
aluminum  under  both  conditions.  In  the  stainless  steel  specimens,  the  measured 
helium  profile  as  shown  in  Fig.  11  agrees  reasonably  well  with  the  predicted  peak. 
However,  the  deuteron  profiles  of  the  two  stainless  steel  samples  are  different. 
The  peak  concentration  of  the  deuteron  profile  measured  in  case  of  a  0.6  mil 
thick  specimen  previously  loaded  with  helium  agrees  well  with  the  predicted 
damage  implant  as  seen  in  Fig.  13.  The  profile  of  deuterium  in  the  sample  not 
loaded  with  helium  was  not  very  statistically  significant  above  background 
effects  as  shown  in  Fig.  14. 

Checks  of  solid  angles,  alignment  of  the  apparatus  and  depth  resolutions 
were  made  with  standard  helium  implanted  samples,  both  surfaces  implants  and 
mid-sample  implants.  Similar  checks  were  made  for  the  deuterium  profiles  by 
use  of  a  CD?  target  5ym  thick  and  with  similar  CD2  specimens  placed  in  front 
and  behind  the  aluminum  and  stainless  steel  samples.  The  surface  specimens 
provide  data  to  correct  for  multiple  scattering  effects  occurring  after  the 
primary  elastic  event. 

The  measured  peak  values  of  the  respective  profiles  are  shown  in  Table  3  in 
comparison  with  the  predicted  implanted  values.  While  a  significant  difference 
between  the  behavior  of  aluminum  and  stainless  steel  is  evident,  the  experiment 
will  be  repeated  with  implants  at  least  an  order  of  magnitude  higher  dosages. 


TABLE  3 


PEAK  CONCENTRATIONS  OF  IMPLANTED  AND  MEASURED  PROFILES 


Material 

Thickness 

(Mils) 

Implanted  (appm) 

He  d 

Measured  (appm) 

He  d 

A1 

1.2 

12 

6 

10+3 

5+1 

A1 

1.2 

0 

7 

- 

6. 5+2. 5 

ss 

0.6 

22 

10 

20+5 

10+3 

ss 

0.6 

0 

12 

- 

<4 
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Fig.  13 

Deuterium  profile  for  the  0.6  mil  stainless  steel  sample  implanted 
with  both  helium  and  deuterium  with  the  predicted  curve. 
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CONCENTRATION  OF  DEUTERIUM  (appm) 


Fig.  IU 

Results  of  a  deuterium  profile  of  the  0.6  mil  stainless  steel  sample 
implanted  with  only  deuterium  compared  with  the  predicted 
curve. 


B.  X-ray  Diffraction  Measurements  of  Low-Level  Radiation  Damage 


(L.  S.  Birks,  J.  W.  Sandelin,  and  D.  B.  Brown 
X-ray  Optics  Branch,  Radiation  Technology  Division) 


Background 


1.  It  is  well  documented  in  the  literature  that  x-ray  diffraction  can 
be  used  to  measure  dislocation  density  in  crystals  and  therefore  should  be 
useful  in  measuring  early  stages  of  radiation  damage  in  reactor  materials.  We 
proposed  that  alloy  compositions  which  are  resistant  to  void  swelling  might 
show  different  diffraction  characteristics  than  alloys  which  are  sensitive  to 
void  swelling  (viz,  lower  dislocation  density  because  of  better  recombination 
of  interstitials  and  vacancies).  If  such  a  correlation  could  be  established,  then 
x-ray  diffraction  would  be  an  effective  rapid  screening  method  for  testing 
proposed  alloys  because  the  irradiation  time  would  be  reduced  by  orders  of 
magnitude. 

For  the  initial  experiments  we  planned  to  use  Cu  crystals  with  selected 
low-concentration  additives  whose  specific  effects  on  void  swelling  have 
already  been  reported.  Single-crystal  rather  than  polycrystalline  material  is 
preferred  because  it  allows  more  accurate  x-ray  measurements  in  a  double¬ 
crystal  spectrometer.  The  crystals  were  to  be  prepared  by  Code  5220. 

Progress 

Results  to  date  have  been  limited  to  a  crystal  of  pure  Cu  because  it  is 
the  only  specimen  which  has  been  prepared.  Figure  15a  shows  the  rocking  curve 
obtained  for  that  specimen  as  received.  The  area  under  the  rocking  curve 
(called  the  integral  reflation  coefficient,  R)  corresponds  to  an  initial  disloca¬ 
tion  density  of  <10  /cm  ,  Table  4.  After  irradiation  by  2  x  10~z  coulombs  of 
3.6  MeV  alpha  particles,  Fig.  15b,  the  central  portion  of  the  rocking  curve 
seems  unaffected  but  the  base  is  broadened  considerably.  We  interpreted  this 
as  an  indication  that  the  damaged  layer  did  not  extend  throughout  the  region 
contributing  to  the  x-ray  diffraction  and  that  the  curve  of  Fig.  15b  represents  a 
combination  of  the  narrow  rocking  curve  from  the  undamaged  material  plus  a 
broad  rocking  curve  from  the  thin  damaged  layer.  Therefore,  a  different 
portion  of  the  crystal  was  irradiated  with  a  graded  series  of  al^ha-particle 
energies  ranging  from  1.2  to  4  MeV  and  a  total  fluence  of  2  x  10~z  coulombs. 
This  resulted  in  the  curve  of  Fig.  15c  which  is  broadened  compared  to  Fig.  15a 
but  re^ains2a  nearly  Gaaussian  shape.  The  dislocation  density  has  increased  to 
3  x  lO'/cm  ,  Table  4.  Ne^t,  a  second  irradiation  was  added  to  the  same  area 
as  for  Fig.  15c  (2  x  10"z  coulombs,  graded  energy).  This  increased  the 
dislocation  density,  Table  4,  and  resulted  in  the  rocking  curve  of  Fig.  15d. 

When  dislocation  density  is  plotted  vs  irradiation  in  Fig.  16  it  appears 
that  the  damage  increases  at  a  faster  than  linear  rate  with  irradiation  initially 
but  that  it  is  beginning  to  approach  a  linear  rate  after  0.04  C/cm  . 
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INTENSITY 


DIFFRACTION  ANGLE 
c  (sec.  of  arc) 

Fig.  15 

Rocking  curves  for  Cu  single  crystal.  Irradiation  with  alpha 
particles  of  a  single  energy  (b)  broadens  the  base  of  the  curve  but 
alters  the  central  peak  only  slightly.  Irradiation  with  alpha 
particles  of  graded  energies  broadens  the  whole  rocking  curve  (c 
and  d). 
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TABLE  4 


DIFFRACTION  CHARACTERISTICS  FOR  COPPER  (Cu) 


Run 

Conditions 

Alpha 

Irradiation  _ 

MeV  Coulombs/crn 

Integral 
Reflection 
Coefficient, R 

Dislocation 

Density 

a 

b 

As  Grown 
Irradiated 

3.6  2  x  10"2 

1.03  x  10"’ 
2.65  x  10"’ 

< 106/crn2 
2.5  x  107 

c 

New  Area 

1st  Irradiation 

Graded  2  x  10~2 

1.2  to  4 

2.74  x  10"’ 

2.5  x  107 

d 

Same  Area  as  c, 
2nd  Irradiation 

1.2  to  4  2  x  10"2 

5.17  x  10"’ 

3  x  108 

Conclusions 

Results  of  the  Cu  crystal  show  that  low  levels  of  irradiation  are  easily 
observed  by  x-ray  diffraction  and  that  increments  of  radiation  damage  are 
measurable  quantitatively.  The  limits  of  dislocation  density  which  cap  be 
me^ure^  with  the  x-ray  apparatus  available  range  from  about  10b  to 
101U /cm  .  Thus  the  method  appears  to  have  sufficient  latitude  to  distinguish 
the  effects  of  alloying  additions  to  the  Cu. 

On  the  discouraging  side  of  the  picture  it  appears  that  difficulties  in 
production  of  single  crystals  precludes  carrying  out  the  experiments  as 
originally  planned.  Therefore,  if  the  project  is  to  continue,  we  must  find  a  way 
to  utilize  polycrystalline  material  or  go  outside  NRL  for  the  specimens. 
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